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Abstract 

We analyse statistics of the real eigenvalues of gl(N, R)-valued Brownian 
motion (the 'Ginibre evolution') in the limit of large N. In particular, we 
calculate the limiting two-dimensional correlation function of spin variables 
associated with real eigenvalues of Ginibre evolutions. 



1 Introduction 



Let gl(N, R) be a linear space consisting of all N x N real matrices and equipped 
with the Euclidean inner product 

(A,B) = tr(AB T ). 

The paper is dedicated to the study of gl(N, R)-Brownian motion (M t : t > 0), a 
Gaussian process with values in gl(N, R) with zero mean and covariance 

E ((A, M t ) (B, M a )) = ^(A, B) (a At), A, Be gl(N, R). 
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In other words, (M t : t > 0) is the process of iV x N matrices whose entries are 
independent one- dimensional Brownian motions. The one-dimensional density of 
gl(N, R) Brownian motion with respect to Lebesgue measure is Gaussian, 

j t (A) = (7rt)-£ C -*<^, 

and, for t = 1, this defines the probability measure for real Ginibre matrix en- 
semble [15]. We will therefore refer to the process (M t ,t > 0) as the real Ginibre 
evolution, or simply the Ginibre evolution. As was shown in |22j, the one di- 
mensional distribution of real eigenvalues for the Ginibre evolution converges as 
N — > oo to the one dimensional distribution of particles for annihilating Brow- 
nian motions on R, under a suitable entrance law. In order to understand this 
unexpected connection deeper, one needs to study the Ginibre evolution beyond 
the one-dimensional distribution. As it turns out, such a study can be simplified 
if carried out in terms of spin variables associated with real eigenvalues which we 
will define as follows. For an JVxiV matrix A let A A be a measure on R: 

A A (a, b) = Number of real eigenvalues of A lying in (a, b). 

The spin variable associated with A is a function s(A) : R — > {±1}: 

s x (A) = (_i)A A (-°^), x G R. (1) 

Note an analogy between the spin variables ([1]) and spins in a one-dimensional 
spin chain with real eigenvalues playing the role of domain walls. The following 
elementary remark provides a tool for computing correlation functions of spin 
variables: the spectrum of a real N x N matrix A consists of real eigenvalues and 
pairs of conjugated complex eigenvalues. Therefore, it is easy to verify, when no 
eigenvalue has real part x, that 

s x (A) = sgn (det (A - xl)) . 

When computed with the help of Householder transformations, the correlation 
functions of spin variables reduce to the correlation functions of characteristic 
polynomials for the Ginibre evolution. The latter problem has been thoroughly 
investigated, at least in the context of real Ginibre matrix ensembles [I], |19j . 
All multi-time densities for real eigenvalues can be restored from the multi-time 
correlation functions of spin variables. Namely we have the following relation: 



Ptl,t-2,...,t„\ X ll X 2 , ■ ■ ■ , X n ) 

(\\ n f n d \ ( n 
~2j (n^) E ( X\ s *m{M tm )s Xm+ym {M tn 



(2) 

ym=+0, m=l,2...,n 



where pt 1 ,t 2 ,...,t n (xi,X2, ■ ■ ■ ,x n ) is the n-dimensional probability density for finding 
a real eigenvalue in each of the disjoint intervals dxk at time k = 1,2, ...n 
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(see [22] for details of fl2]) in a special case). From a purely technical point of view, 
it is also useful to consider modified densities defined as follows: 



Pt 1 ,t 2 ,...,t n \ X ll X 2, ■ ■ ■ ,X. 



n 



)l[dx k =E[ l[s Xk {M tk )A M ^{dx k ) 



k=l \k=l 




(3) 



Spin correlation functions can be restored from the modified densities by n-dimensional 
integration (see (1641)). 

In this paper we use the machinery described above to calculate exactly the 
two-dimensional correlation function of spin variables, namely ~E(s x (M t )sy(M t j rT )) . 
This is done by a lengthy exact calculation and is an admittedly modest step to- 
wards the complete understanding of the Ginibre evolution. However we hope that 
some of the structure will be useful in calculating higher dimensional correlation 
functions, and we illustrate this by re-examining multispin correlations at a fixed 
time. 

The paper is organised as follows. In section [2] we state and discuss our main 
results: the large N limit of the two-dimensional correlation function and a new 
integral representation for one-dimensional multi-spin correlation functions of real 
eigenvalues. In section [3] we present the main steps of the computation of the two- 
dimensional correlation function for fixed N, based on the Householder transform 
and the technique of integrals over anti-commuting variables (Berezin integrals). 
In section S] we examine the N — > oo limit. In section we test the techniques 
developed in section [3] to re-derive the result of [2] for the large N limit of fixed 
time correlation functions pi,...,i(a;i, . . . , xk) and pi,... t i(xi, . . . , xk), that is the ex- 
pressions and (|3J) where all times ti, . . . , t# equal 1. 



2 Discussion of results 

Our first result is the de-correlation of spin variables in the real Ginibre evolution 
for a fixed positive time lag. 

Theorem 1. (Propagation of temporal chaos.) For x, y e R and t, r > 



This shows the large- N and small time lag limits do not commute. It suggests 
that in the large- N limit, the system of real eigenvalues in the real Ginibre becomes 
memoryless, that is marginals at distinct times become independent. This result 
supports the belief that real Ginibre ensemble belongs to the class of systems with 



lim E(s x (M t+T )s y (M t )) = 0. 
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long range interaction. For a long-range first order system of N particles, the 
scale of temporal correlations is expected to be of the order of 1/N, see [3] for 
a good heuristic discussion of the Hamiltonian case. To study these short-scale 
correlations we consider the scaling limit r = j£ and N — > oo. Therefore our 
second result concerns the behaviour of spin-spin correlation function in the real 
Ginibre evolution for time lags of order 1/N. 

Theorem 2. (Decay of temporal correlations.) For x, y G R and t, r > 0. 

Jim E(s x (M t+j _)s y (M t )) = 
where erfc is the complementary error function. 




The large- N limit of the two-time spin-spin correlation function characterized 
by Theorem [TJ and Theorem [2] is different from the two-time spin-spin correlation 
function for the system of one-dimensional annihilating Brownian motions obtained 
in [21] . Unfortunately, this means that the conjecture concerning the law of real 
eigenvalues in real Ginibre ensemble made in [22J is false. 

Our final result concerns a novel integral representation for the modified K- 
particle density of real eigenvalues at a fixed time. 

Theorem 3. ( Ginibre ensemble and anti- self dual Gaussian symplectic ensembles.) 
Fix X\ < x 2 < . . . < xk £ R- Then 

lim pi,...,x(xi, x 2 , . . . , x K ) = CVA(x) / fi H (dU)e-^ Tr ( H - HR )\ 

N ~>°° ' JU(K) 

where Ck is a positive constant, H = U XU< is a Hermitian matrix with eigenval- 
ues Xi,X2, ■ ■ ■ ,xk, lift is Haar measure on the unitary group U{K), H R = JH T J 
is a symplectic involution of matrix H , J is the canonical symplectic matrix, and 
A( x ) = rii>j( x « ~ x j) i s the Van-der-Monde determinant. 

In the present paper we only use Theorem |3] to show how the formalism of spin 
variables can be employed to re-derive the Pfaffian expressions for the correlation 
functions of real eigenvalues in the real Ginibre ensemble. These were originally 
obtained in [2]. However we suspect that there is a generalization of the Theorem 
|3] to the case of multi-time correlation functions which should be useful for the 
study of multi-time statistics of the Ginibre evolution. 
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3 Calculation of the two-dimensional spin-spin 
correlation function. 



The calculation detailed below draws on the ideas of Edelman et. al. |6J and 
Sommers et. al. [20] to compute the two-dimensional spin-spin correlation func- 
tion while avoiding the mathematically difficult transition to the eigenvalue rep- 
resentation. Edelman used Householder transform to compute the density of the 
eigenvalues in real Ginibre ensemble while Sommers employed Berezin integrals to 
simplify his derivation of the kernel of the Pfaffian point process associated with 
the real Ginibre ensemble. The key fortuitous link which makes our calculation 
possible is that s x (A) = sgn(det(/i — xl)) (see eq. ([1])) which combines nicely with 
the Jacobian | det(A — xl)\ of the Edelman transformation (see below) to produce 
the characteristic polynomial det(A — xl). The latter is easier to average over the 
Ginibre evolution. 

Instead of computing the spin-spin correlation function directly, we compute the 
time-space point modified density function first: 

p N (t,y;t + T,x)dxdy = E [s y (M t )A Mt (dy)s x (M t+T )A Mt ^(dx)] . (4) 

The notation stands for densities where we make sure here to take the limiting 
infinitesimal intervals dx and dy to lie just to the right of the points x and y 
respectively. Note that the function x — > E[s a (M t )A Mt (db)} has a discontinuity 
as a crosses b, indeed the density flips sign. Given p^, the spin-spin correlation 
function can be restored by integration: 

/x r-y 
dx' dy' p N (t,y';t + T,x'). 
-oo J — oo 

It follows from the definition (T4]) of the modified density that 
Piv(t, y, t + r, x)dxdy 

= [ dM [ dHs y {H)A H {dy)s x {dM)A M {dx) lt {H) lT {M - H). (5) 
Jr n2 Jr n2 

The plan of attack is, via changes of variable that exploit symmetry of the inte- 
grand, to reduce the dimension of this integral. Indeed by the end of this section 
the integral over R 2Ar is reduced to a one dimensional integral over [0, 1]. 

The representation of dM integral in Edelman coordinates. The trans- 
formation used by Edelman et. al. in [6] is for an Af x iV matrix M with a real 
eigenvalue x and corresponding eigenvector v G the upper half of the N — 1 

dimensional unit sphere in H N : 

M = P,.M' I\. 
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where P v is the Householder transformation [TT] that reflects in the hyperplane at 
right angles to the vector v — (where is the unit vector (0, . . . , 0, 1)), and 
M e is in Edelman block form 



M e 

T 

W X 



with Mq an (N — 1) x (N — 1) real matrix, w G R 1 and x G R. Later we need 
the explicit form 

Pv = I _^-^~^) T (6) 
1 - v N 

where vn is the iVth co-ordinate of v. Note that P v is orthogonal and symmetric, 
P- 1 = PJ = P v . Let 

F(M)dy= [ dHs y (H)A H (dy) lt (H) lT (M - H). 
Jr n2 

Fix t, r > and x,y throughout and write as shorthand for p N (t,y;t + r,x). 
Using the known expression for the Jacobian of the Edelman transformation [6], 
we rewrite the expression (jSJ) for p N as 

Pn = I dv [ dMl [ dw\det(M^ - xI)\s x (M^F(P v M e P v ) 
= [ dv f dMl I dwdet(M% -xI)F(P v M e Pv) 

Js+r^ JR(N-i) 2 JrN-i 



where dv stands for surface measure on the sphere. We have also used s x (P v M e P v ) = 
s x (M e ) since s x (A) depends only on the eigenvalues of A, and s x (M e ) = s x (Mq) 
since we took limiting densities where dx just to the right of x. Note the disappear- 
ance of the absolute value sign in the last expression leaving us the more tractable 
determinant. 

Evaluation of the dv integral. We claim that F(0 T MO) = F(M) for an 
orthogonal matrix O. To see the latter note that the Gaussian density has this 
invariance and therefore 

F(0 T MO)dy = [ dHs y (H)A H (dy)j t (H)j T (0 T MO - H) 

dHs y (H)A H (dy) lt (H) lT (M - OHO T ) 

[ dHs y {0 T HO)A° THO {dy) lt {H) lT {M - H) 
Jrn 2 

F(M)dy. 



\ - 
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Here we used that s y (H) and A H (dy) depend only on the eigenvalues of H and are 
therefore invariant with respect to the transformation H — > T HO. Then 

Pn = I dv f dMl [ dwdet(M%-xI)F(M e ) 

JS'+_ 1 JR(N-1) 2 Jrn-i 

= hs N -i\ [ dMl [ dwdet(M* -xI)F{M e ) (7) 

2 Jr(n-i) 2 JrN-i 

since the integrand is independent of v. Here |Sat-i| stands for the surface area of 
the unit (N — 1) dimensional sphere living in R . 

Evaluation of the dw integral. To integrate over w we must express F(M e ) in 
terms of Mq and w. Let us represent H in block form: for z G R and a, (3 e R^ -1 , 
write 

H=( H }I> 
y a z 

Expanding the Gaussian densities we find that F(M e )dy is given by 



fvr 2 ' 



Hr)-^ [ dHs y (H)A H (dy) e -(H^^) e -±tr(M*M-) e ftr(M^)_ (g) 

The traces expand in block form to 

tr(M e M eT ) = tr(M e M eT ) + |w| 2 + x 2 , 
tr(M e H T ) = tr(M e # T ) + w T a + xz. (10) 

Substituting (jUj) and f lTUj) into (j7j) we find the following representation for p; 

p N dy = i|Sjv_i|(7r 2 tr)~^e~^ / dM e [ dw [ dH 

s y (H)A H (dy) det(M e - rrl) e ~(^) tr( ^ T) 
e -i(h[2-2^ a _ 2M ) e -itr(M«M ^) e ftr(M«HT)_ 

The dw integral is then computed using a standard formula for Gaussian integrals: 

d«-^H ! - 2A ) = H^e^l ! . (12) 

R N-1 

Evaluation of the dM$ integral. To compute the dM$ integral we first re- 
express the determinant as a Berezin integral over anti-commuting (Grassmann) 
variables: 

det(M e -xl)= [ d<pd~4> efm-xi)*. 

J R 0\2(N-1) 

This leaves the following dM^ integral in f TTTj) 

h := / dMl e -^mMf) e ftr ( M ^)) ^m*^ (13) 

JR(N-1) 2 



This is a Gaussian integral with identity covariance matrix but with the linear 
term of the exponent depending on anti-commuting variables. Nevertheless the 
rules for this integral are as if it were a standard Gaussian (see [32] for details), 
and the value can be found by locating the critical value of the quadratic form in 
the variables M^(i, j). The integrand in (fl~3l) is the exponential of 

- itr(M e M eT ) + -tr(M T T ) (14) 

T T 

where T = H + ^(fxf) T . The critical point of this quadratic is = T. Substituting 



this back into (fT4j) one obtains, noting that 



T\T 



,-rT 



+ I tr( rr T ) = + l -tr{{H,+ T -<p<f){HT- W)) 

T Til 

= +-ti(H H^) + ti(H W) 

T 

This gives the value 

h = (7rr )^^ tr ^oW/W (15) 

Substituting in the dw and dM^ integrals ( !T2|) and ( IT5|) into ( TTTj) we reach 

1 N x' 

PNdy = -|S , Ar_i|(7rt) _ ~(7rr) _ ^"e _ ~ / dH 



' R n2 J_R0|2(JV-1) 

s y (H)A H (dy) e ^^ + ^ tT{HHT) e H 2xz +\ a \ 2 ) e ^i{H H^)+t '(H -xi)<f> _ 

The integral at hand can be simplified by an orthogonal transformation on the H 
variables. To simplify the implementation of the transformation, we will re-write 
the above expression in terms of matrix H rather than its sub-matrix Hq. To this 
end, we use an extended set of iV-dimensional Grassmann variables ip, ip that agree 
with in the first (N — 1) co-ordinates, 

ip = {0,?M, V> = {4>^n}- 

Then 

tp T (H - xl)ip = (f) T (H - xl)4> + (4> T /3)ip N + ifj N (a T (f)) + (z- x)tfj N ip N , 
and we can re-write 

/ d(f)d0e* T{Ho - xI),t> = [ di)d^^ N ^ N e^ T{H - xI) ^ 

J R 0\2{N-1) JR0\2N 

since the integral over the pair dipNdip N can be done first on the right hand side 
and the term ip N ipN forces it to take value one. Using also ti(HH T ) = It(HqHq) + 
|/3| 2 + \a\ 2 + z 2 this leaves 

~p N dy = -\S N ^\(7rt)--(nTy-e-- / dH s y {H)A H \dy) 



N 2 



##^iv e^*-*™ e -^ HH ^ e -m 2 ^ 2 -^). (1 6 ) 

R 0\2N 



Representation of the dH integral in Edelman variables. A second Edelman 
change of variable is H = P v H e P v where 



H e 

T 

w y 



We need to reconcile this with our earlier representation flSJ) of H in block form. 
A lengthy but straightforward calculation shows 

z = (v T H^v) + v N (w T v) + yv 2 N , 
(3 = HqV — (v t Hqv)(1 — vn^v + (w T v)v + yv N v, 
\P\ 2 + z 2 = (v N y + w T v) 2 + v T HfH*v. (17) 

Next 

^ T (H - xl)ip = (P v ^) T (H e - xI)P v ip, 

which suggests a change of Grassmann variables <fi = P v ip, = P v ip- Under this 
change we find with the help of (EJ) that ip^ = v T 4> an d V'jv = v T (j). Therefore, 
Berezin integral in the right hand side of ( 1161) transforms as follows: 



T, rr „ , / „ —rT , 



I 2 := dipdip ^ N ip N ^ {H ~ xI)lp = / dcPdcP^^^^e 4, {He ~ xI) t (18) 

J R 0\2N JR0\2N 

Using this and the Edelman substitution in (|T6|) we reach 

Pn = -|Sjv-i|(7r£) 2 (ttt) 2 e T I dv dH^ / dw 



J N-1 

-T , 



>R(N-i) 2 Jrn-i 

d(f)d4> (y T 4>)(v T <P) ( He -*^ det(# e - yl) 



R Q\2N 

e -\tl{H-H- T ) e -±((v N y+w T v) 2 +v T HfH$v) e ^(v T HSv+v N (w T v)+yv%) _ 

To prepare for the integration over w and Hq we found it convenient to integrate 

out some Grassmann variables. First we let <j), <p be the restriction of 0, <fi to the 
first (N — 1) co-ordinates so that 

f{H e - xl)(j> = JVo - + l> N {w T 4>) + (y- x]4> N ct> N . 

We will integrate out the final pair of variables <pN, <Pn- Extracting the terms 
depending on 0tv, 4>n gives 

N d$ N (v T J + v N N )(v T 4> + v N <f> N )e*»^ T +') + <*- x 'ft"* N 

R 0\2 

d(f) N d(f> N (v T $ + v N ~(]) N )(v T (f) + v N (p N )(l + lf> N {w T <j>)){l + (y- x)<j> N <j> N ) 

R 0\2 

v 2 N + (y- x){v T $){v T 4>) + v N {w T $){y T J). 
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This leaves us just with the integrand depending on shortened anti-commuting 

variables 0, <p and so we can drop the hats in the notation for integration variables 
replacing them with (f>, (f> e R !^ -1 ). Thus we have shown that the Berezin integral 
(|T8|) has become 

d<pd~4> (v 2 N + (y- x)(v T 0)(v T <j)) + v N (w T (j))(v T 0)) e^^o-^ 

R 0\2(N-1) 

(v% + (y- x)({) T 0)(£ T 0)) e^ (wT</,)( ' T?) f^-^t 

i?0|2(iV-l) 

Using this, representing det(ifQ — yl) as a Gaussian integral over a second set of 
Grassmann variables ip, ip an d expanding 

tr(H e H eT ) = tr(H^Hf) + \w\ 2 + y 2 

we reach the following expression for p^: 

1 _n2 _N f f f 

Pn = -|Sjv-i|(7r£) 2 (ttt) 2 e T I dv I dH$ I dw 

d^d^ e r(HS-xi)<P+r(HS-yi)4> 

RPMN-1) 

+ (y _ x )(v T Mv T <P)) e ^ {wT ^ T ^ e -\^(H^ )+ \ w fW) 

e -±((v N y+w T vf+v T HS T HSv) e ^(v T HSv+v N (w T v)+yv 2 N ) _ ng\ 



Evaluation of the second dw integral. Collecting all the terms in the integrand 
of f lT9|) depending w one reaches the second dw integral, which happens to be 
Gaussian: 

j 3 = I dwe~^ H2 e -^ 2 (y- x ) v N{w T v)+{w T v) 2 ) e ^(w T <p)(v T ^)^ ^ 

The exponential has a quadratic term t~ 1 w T Mw with covariance matrix 

M = I+-vv T . (21) 

T 

Note that 

M' 1 = I t ——vv T = I-—vv T where D := det(M) = 1 + -\v\ 2 . (22) 

T + t\v\ 2 TD K J 

The integrand in (I20j) is the exponential of 

1 2 

- -w T Mw + -w T 7 (23) 
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where 7 = ^(x — y)v^v — ^^{v 7 '"</>) <f>- The critical value of w in then M lr y and 
substituting into ( 1231) one obtains t~ lr y T M~ lr y which is 

~{-{x - y)v N v T - -L(v T ^ T )M-\-(x - y)v N v - -^({) T 0)0) 

t T ZVn T ZVn 

= -^{x- yfv 2 N v T M~ x v - ~(x - y)(v T '(j))v T M' 1 (j). 

Using the form ( l22|) for M _1 one finds that v T M~ 1 = D^ 1 v T and this leads to the 
answer 

/ 3 = £)-5(7rt) V e -^(-J/)(^)(- T ^) e^-f) 2 ^! 2 . (24) 

Evaluation of the integral. Using expression fl2TT) for matrix M we can 
verify the formula 

-MH^Hf) + -v T HfH e v = -MH^MHf). 

t T t 

Combining all the ifg-dependent terms in the integrand of (Tl9|) and using the above 
identity we reach the integral 

h= f dH e e ^oW T ^ e -^(H S MHf) e ^ H§ ^ (25) 

The quadratic part of the exponential tr(HQMHQ T ) can be considered as a quadratic 
form in the variables (iJg (*? j) : hj — 1, ■ ■ ■ , N — 1). The corresponding covari- 
ance matrix has determinant Z)^ 1 . The easiest way to see this is to take v = 
(1, 0, . . . , 0) and examine the effect of the transformation Hq{i,j) — > (HqM)^, j) 
for each i, j. The integrand in ff25l) is the exponential of 

- hr(H e MHf) + - t tr(H*r T ) (26) 



where 



r = + + —vv T . 

I I T 



The critical value of Hq is = TM 1 . Substituting into ( 12 6 p one obtains 
t~ 1 tr(rM- 1 r T ) which is 



— ti(vv T M- l vv T ) - -tr( (f)(f) T M~ 1 ^ ) + — tr( 6d> T M~ 1 vv T ) + — 

T 2 T T 



and evaluates, using the form (|22|) for M 1 , to 
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This gives 

J 4 = D -(N-l)/2^(N-lf/2 e -^ T ^ T ^) e ^(fi T M* T V)(^?) 

e &((« T ^(fi T 0+(^)(« T i0) e ^H 4 . (27) 
Substituting ([MD and ([27]) into (PSD we reach 

-ISiV-iK^r)-^^ /" ^ jD -f e -^( a: 2 -2^ :r2/ +^(i +T t-i)) 

— ' J 'S'at t 



-tT, -t-T, 



d<j)d<j)dil)dip (v 2 N + (y- x)(y T (f)){y T (f))) e~** * 

_R0|4(AT-1) 

Integration over Grassmann variables. Let R v : R^ -1 — >■ 1 be an or- 
thogonal transformation that sends v to (\v\, 0, . . . , 0). We also change variables 
<p,ip,<p,ip to R v (j), R v tp, R v (j), R v ip. Under this change i) T (j> becomes \v\4>i e.t.c. In 
terms of new variables the integral in (1281) is 

~ PN = -ISn-xUttHt)-^ I WTe-aj^-HWa+rt- 1 )) j 5) (29 ) 

where 

h — dcpdipdipdip (y 2 N + (y — x)\v\ 2 ^(pij 

e ~xlj> T <f>-yil> T ip e -t0 F -4>)((f> T ili) e ^^-(4> T ip)<t>i-ipi e ^-(xip 1 i(i 1 +ylf; i (f> 1 ) 

= v 2 N [ #### i' ' ^ e" x ^-^ 

e -|(0 T V')(0 T V') e %7S"(^ T V')</'i'/'i e ^-(a^i^l+3/0i<frO 

Notice the appearance of fourth-order terms (0 ip){(f) T 'ip) and (0 VO^iV'i in the 



above integral. To deal with these we follow (20] and apply Hubbard- Stratonovich 
transformation to to convert the Berezin integral I to the integral over both com- 
muting and anti-commuting variables but which is Gaussian with respect to the 
Grassmann variables, by using the identity 

e ab = [ ^! e - z W J (30) 

In this identity, we have z = x + iy, z = x — iy, and dzdz is shorthand for a 
Lebesgue integral over R 2 of this complex integrand. The result is 

h = v% [ ^ [ [ 



_ T rp +i"|2 — (y~ x )l^l 2 T ± 
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where a + /3 = 1, 7 + 5 = 1. The final answer should not depend on a, (3,^,5, 
which provides us with a consistency check for the calculations below. The inter- 
nal Berezin integral breaks into the product of iV — 1 independent 4- dimensional 
integrals, of which (N — 2) integrals are identical: 



dzdz 



dwdw 



N-2 t 



where 



and 



/R2 7T J R 2 tt 

#### e (*> a ***-U)^ e &V"M e -***-jH> 



R 0\A 



R 0\4 



Each Berezin integral is exactly in Pfaflian form (that is Pf(-A) = f d£ exp(— £ T A£/2)). 
Hence I e has value 



/0 x -(I)' 



Pf 







\ 








(ir-+(S)' 



V 



y 
o 



/ f 2|0|2\7 



2tD 



wz 



and J7 has value 
/ 



Pf 



x — y '-' ^-LJ- 



\v\ 2 t_±r 
N 



w 







v • 



x — y\v\ 



jt + T 

tD 






y — x \v\ 



y- 



*L21: 



X 



t 2 |t>| 2 



WW 



a ft 2 \v\ 2 \ 5 



wz + I - 



^1^)12X7 



2rL> 



Integration over Hubbard Stratonovich variables. To compute the integrals 
over complex variables z, w notice that 



dzdz f dwdw_ e _^_^ zk ^ wm _ fl 



'zdz r dwo 
^ Jc ^ 
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is zero for integers k, /, m, n unless k = I and m = n. Note that I 7 contains terms 
in Wz and wz. So if we expand (I 6 ) N ~ 2 as 

/ t \ N ~ 2 ( t \ N ~ 3 ftV ft 2 \v\ 2 V 

[xy + -zz) + (N — 2) [xy + -zz) ^ ^±Lj W z + ... 

we may ignore the other terms as they give zero contribution upon integration 
with respect to w,z. Substituting the values of I e and Ij we find 

2 f dzdz f dwdw _ wW 

h ~ VN Jc^rJc^r 

(4 (* " 2/I " |2 t^) ( y ~~ ^ |2 ^) + r~ z ~ ^ wW ) { xy + H 

These integrals can be computed using the formulae: 

- / dzdze~"(a+ \z\ 2 ) N = [ dre- r (a + r) N 
71 Jb? Jo 

=ib(T) ak r e ~ rrN ~ k=EN{a) ( 3i ) 

where En(x) is the scaled and truncated exponential polynomial, 

N k 

Similarly 

- [ dzdz\z\ 2 e~ z *(a + \z\ 2 ) N = E N+1 (a) - aE N (a) 

= (N + l)E N (a) - NaE N -!(a). (32) 

So, integrating over w,w and pulling powers of t out, 

(1 ( rfl t + T \( 112M t - t 2 \v\ 2 \(2xy \ N - 2 
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Finally, integrating over z,z, we find, using (131]) and ( 1321) . 



7V-2 



,t + T 



x — y\v\ 



N-2 .2|-|2„,2 



t z \v\ V 



N 



+v 2 N 



2tD 

N-l 



-D 



En -2 



y — x \ v \ 

2xy s 



-D 



En -2 



2xy 



(N - 1)E N -, 



2xy 



- {N-2) 2 ^E N _ 1 



-v N (N-2)[- 



N-3 ,3 



t a \v\ z 
AtD 



E 



N-2 



2xy 



t 

2xy 



N-2 



x — y\v 
(N-l)v% ft 



D 



jt + T 

tD 

N-l 

E 



y — x\v 
2xy 



N-2 



t 



(N - 2)v 2 N ft 
D 



N-l 



2xy 



En-3 



tD 



f 2fy 
V * 



E 



N-2 



t 

2xy 



E 



N-3 



2xy 
2xy 



Substituting this into (129]) . and using tD = r + t\v\ 2 , gives us the following exact 
representation for the modified density pN in the form of the integral over the 
(N — l)-dimensional sphere 



p N (t,y;t + r, x) 
2 Sjv-i 



t 2 (2vr; 



N 



dv (rr 1 + \v\ 2 )' 



t + r 



x — y\v\ 



T + t\v\ 



t 



T + t\v\ 



N-2 



(N - l)trv 



N 



2(r + t\v\ 2 ) 



E 



N-2 



y — x\v 

2xy\ (N - 2)trv 2 N 2xy 



E 



2xy 



t 



2(r + t\v\ 



t 



E 



N-3 



2xy 



(33) 



1 — v N and so the integral can 



The integrand depends only on % and \v\ 2 
be further reduced, as we do in the following section, to a single integral over 
v n € [0, 1]. Moreover we will obtain very much simpler expressions for both pN 
and the spin-spin correlation function valid in the limit N — > oo. 



4 Analysis of the large N limit of the correlation 
function. 

A direct evaluation of ( 133]) leads to the following conclusion: 

Hm p N (x,t + r;y,t) = 0, 

7V-»oo 



for any r > 0. To study short-scale correlations, we consider the scaling limit 
r = and N — > oo. We will calculate 

p(x,T;y } t) = lim p N (x,t + (T/N);y,t). 

It is possible to calculate p(x, T; y, t) directly. The calculation is somewhat lengthy 
but the answer turns out to be, as expected, translational invariant. Assuming this 
translational invariance, it is possible to calculate p(x,T;y,t) as p(0,T;y — x,t), 
which is a similar but less messy task, and this is the calculation which we present 
here. 



Firstly, when y = the expression (1551) simplifies considerably. Furthermore, 
using v% + \v\ 2 = 1, the integrand is invariant under rotations in the equatorial 
plane and we may exploit the identity 

dvg{y N ) = \S N - 2 \ / dzg{z){\ - z 2 )^ . 



We find 

p N (t,0;t + (T/N),x) = d(N) / dz (1 - z 2 )~ 5/2 e wwi^y 



o 



T V^' 1 / 2 2Nx 2 (l-z 2 " 
1 + Nt(l - z 2 ) J V (N-1)T 



where 

, Ar , (JV-1)!._ „_ . T 
d(JV) = \ n ' \S N ^\\S N „ 2 \ 



Using the volume formula [*Sjv_i | = f£ — . , and asymptotic T(z) ~ z z z x / 2 e z for 

large z, it is straightforward to check that C\(N) — > T/(irt 2 ) as N — > oo. Moreover 
one may justify the passage of the limit inside the dz integral, for example by 
dominated convergence using 

JV 

This limit gives 

T f 1 -, , 9 , = „ r / , 2x 2 (l-z 2 )\ 

p(t, 0; T, x) = — y (1 - 2 2 )- 5/2 e e «(i=^) L_ L . 

The substitution (1 — z 2 ) = (1 + u?) _1 and = |w _1 / 2 (l + w)~ 3 / 2 dw yields 
p(t,0;T,x) 

I -(4+1) r dwe -(4+^ ( w i/2 _ 2 ^ w -i/2 



2nt 2 J \ T 

T e -W + i)flff! + lV 3/2 _^fE! + lV 1/2 | (34) 



20rt 2 V 2 V * 2 V T \ t 2t 
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which is the final answer for the scaling limit of the modified two-point density in 
the large N scaling limit. Note that the limit 



is consistent with the one-dimensional correlation function obtained in [2] (the 
presence of the delta function at zero is expected, see [21] for an explanation for a 
similar two-time density). Let 



p(t,y;T,x) = lim E (s y (M t )s x (M t+ T)) 



be the scaling limit of the unmodified two-time spin-spin correlation function. This 
is given by 



p(t,y;T,x + y) = lim p N (t, y; t + (T/N), x + y) 

AT— »oo 

roc rx+y 



poo r x +y 

= - lim 4 / dy' / dx' p N (t,y';t + (T/N),x') 

N ^°° Jy J-oo 

poo r x ~\-y 
= -4 / dy' / dx' p(t,y';T,x') 

J y J — oo 

= -4 / dy' dx' p(t, 0; T, x' - y') 

J y J — oo 

/x 
dz{z-x)p(t,0;T,z). (35) 
■oo 

Justifying the passage of the limit iV — > oo inside the integral in (135]) needs a little 
care, and we delay the argument until the end of this section. Substituting ([34 
into the above, we find that 



4 r°° f 1 

p{t, y- T, x + y) = -j=- t J dz{x - z)e~ a ^ i^-ra^ 2 (z) - (a(z) - r) a - l ' 2 [ 



z 



where r = a{z) = ^-+r. Since p(t, 0; T, x) is even in x so too is p(t, y; T,x + y). 
Taking x > 0, the change of variables w = a(z), so that x/t 1 ^ 2 = (a(x) — r) 1 / 2 , 
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gives 



p(t,y;T,x + y) 

2 r° , ... / 1 



dwe w ( -rw 3/2 — fty — rlio 1/2 



K J a (x) \2 

vntj a (x) \2 



-I/ 2 + -1= / dwe-^Va 

V 71 Ja(x) 



= -= / dwe^ii)- 1 ' 2 , 

V 71 " Ja(i) 

where the last two equalities follow by integration by parts, and all the boundary 
terms cancel. This final integral is the same as the complementary error function 
giving 



p(t, y;T,x + y) = erfc ( \j j + ^ j . (36) 

For T = 0, ( 1361) reduces to a well-known answer for continuous limit of the single 
time spin-spin correlation function in the Glauber model [8]. As T — > oo, the 
spin-spin correlation function decays exponentially, as expected. 

We now return to complete the justification of ( 13 5 p . We need to pass to the limit 
to show 

POO PX POO PX 

/ dy' / dx'p N (t, y'; t + (T/N), x') ^ dy' dx'p(t, y'; T, x'). (37) 

J y J —oo Jy J — oo 

By y — > —y, x — > —x symmetry we may suppose y > 0. Using the fact that pN is 
a mixed derivative we know 

p N (t,x;r,y)dx = / p N (t,x;r,y)dy = 0. (38) 
r Jn 

The same property holds for the limit p(t,x;T,y). When y > and x > it is 
convenient to use fl38l) to rewrite (1371) as 



oo ^0 py rx 

dy' / dx'p N (t, y'- 1 + (T/N), x') + / dy' / da/p^t, y'; t + (T/N), x') 

y J —oo J — oo Jo 

poo p0 py px 

— > / dy' / dx p(t,y';T, x) + / dy' I dxp(t,y';T,x). (39) 

Jy J —oo J —oo Jo 
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This ensures that the region of integration has only a bounded region where xy > 0, 
namely [0, x] x [0, y}. 

We will replace pN by a further modification p^, where we replace each occur- 
rence of the truncated exponentials En(z) by the un-truncated exponential N\e z . 
Thus we define (compare with (1331) ) 

p N (t,y;t + t,x) 



->N-l\ / , / ,_1 , i~i2\-4 TTTr(a; 2 -2fi^+S/ 2 (l+T*" 1 )) 2£E 



9+ 

|2 ^ + T A ( IN2 * 



r + t\v \ l ) \ t + t|f 

(AT-l)Ur^ (JV-2)!*r^2a^/ 



(40) 



2(r + t|{)| 2 ) 2(r + t\v\ 2 ) t 
This does not change the pointwise limit, and we still have 

T; y, t) = lim p N (x, t + (T/N); y, t). 

iV— >oo 

We need a uniform bound for px(t,y\t + (T/N),x). We write C(t,T, . . .) for a 
constant whose value may change but whose dependency is indicated. The key 
term is the exponential: when xy < we may bound 

1 1 , ,-> (x 2 -2v 2 M xy+y 2 (l+TN- 1 t- 1 )) 2m . * 2 +v 2 

and when > we use 

e - w -i + ,H4 ' 2 -HW(i+Tr'r'))^ = e - ^ jW ( (3! - v )a + (y»- 2ay )riy-it-i) 

(^-«) 2 23-y 

< e T +* e * . 

The remaining terms in the integrand in p^(t,y;t + (T/N),x) can be simply 
bounded by C(N,t,T)(l + x 2 + y 2 )\v\~ 2 , and carrying out the integral over 
yields the bound 

( x 2 +y 2 

p N (t,V,t+ (T/N),x) < C(t,T)(l + x 2 + y 2 ) e ~™ 2 „, iixy -°> (41) 



_{x-y)_ 2xy . 

e T+t e t if xy > 0. 



When y > and x < the limit 



/*00 /*£ /»00 y^X 

/ dy' dx'p N (t,y';t + (T/N),x') ->> / dy' / dx' p(t,y';T,x') 

Jy J —oo Jy J — oo 



follows from ( 14T|) by dominated convergence. When y > and i > we switch to 
the form ( 139]) and obtain, also by dominated convergence using ( 14T1) . 

dx'p N (t, y'\ t + (T/N),x) + dy' dx'p N (t, y'\ t + {T/N), x') 

-oo J — oo Jo 



— >• 



/ dy' / dx p(i, y';T, x) + / dy' / dxp(t,y';T,x'). 

Jy J —oo J —oo Jo 
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It remains to show the error from approximating pjy by pN is negligible. However 
this is a simpler task, using the fact that E N (z)/N\ is close to e z on a ball of radius 
N 1 / 2 , and using the simple bound \En(z)\ < iV!e' 2 ' outside this ball. We omit the 
details. 



5 Fixed t multi-point intensity functions. 



As another application of the formalism of spin variables we will show how to re- 
derive the result of Borodin and Sinclair [2] for all the one-dimensional densities of 
real eigenvalues of the real Ginibre ensemble in the limit N — > oo. Another reason 
for going through this calculation is to prepare ground for the future study of 
multi-time correlation functions. Let us fix time t = 1 and a set of multiple space 
points x\ < X2 < ■ ■ ■ < xk- The answer for an arbitrary time t can be obtained 
from the t = 1 answer by the diffusive re-scaling x — > As we have done before, 
we will first compute a modified density p^ = Pn( x i, x 2i ■ ■ ■ > x k) defined by 

p N (x 1 ,x 2 , ...,x K ) dxi . ..dx K = E N I Y[s Xk (M 1 )A Ml (dx k ) 

\k=l 

where where the subscript iV on the E^v means that the averaging occurs over the 
N x N Ginibre distribution. As before, we choose right hand limits for the density, 
that is where the intervals dx k denote infinitesimal intervals just to the right of 
the point x k - 

The answers are zero for odd K, so we take even K throughout. Moreover 
it is convenient to consider only even iV throughout (which avoids us tracking 
various ± signs). 

Equivalence with a correlation function of characteristic polynomials. 

The integral over the Gaussian density 



K 

p N {Xi,x 2 ,...,x K )dxi...dx K = / dM^i{M)Y\ s Xk (M)A M {dx k 

J R N2 Li 



can be treated using the Edelman transform for the eigenvalue lying in dxx as in 
Section [31 and after integrating over the half sphere S^-_ x , we obtain 



Pn( x i, x 2, ■ ■ ■ , x K ) dxx... dx K -i 



K-l 



hSn^-^e'^n^ (det (M.-xkI) J] s Xk (M 1 )A Ml (dx k 



k=l 
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Another Edelman transform about the eigenvalue lying in dxx-i yields 



Pn{xi,x 2 , x K ) dxi... dx K -2 
1 
4 



1 N—l N — 2 2 2 

'SjV-l SV-2 7T~~ - e - x K- x K-i ( Xk _ x - x K ) 



K-2 

E N ^ 2 ( det (Mi - a*I) det (M a - J] s Xk {M 1 )h M "{dx k ) 



fe=i 



A further (it" — 2) applications of Edelman transform will lead to the following 
expression for the modified density: 



Pn{Xi,x 2 , . . .,x Kj 

K / K 



A(xl 



,- )K J] (iSN-kW-^e-^ E N _ K J] det (Mi - x m J) (42) 

fc=l \m=l / 

where A(x) = Yli<i<j<K( x j ~ Xi ) * s ^ ne Van-der-Monde determinant. Therefore, 
the problem of computing the modified density has been reduced to the computa- 
tion of the expectation of the product of characteristic polynomials of the random 
matrix Mi. 

Integral representation for product of characteristic polynomials. Such 
a computation has been carried out in [2D] by exploiting Berezin integrals. As it 
is rather close to the calculation already detailed in Section [3] we will just present 
the final answer. 

K 



E N ] J det (Mi - x m I) 



V7Tt=l 



n 

i<p<q<K 



dz pq dz pq ^_\ Zpq \2 
R2 7T 



Pf I ?x ^ I • (43) 



Here each dz pq dz pq is shorthand for Lebesgue measure on R 2 and arises from re- 
peated use of the Hubbard Stratonovich transform; the matrix A is a diagonal 
K x K matrix with entries (xi,x 2 , ■ ■ ■ xk)', and Z is a skew symmetric complex 
K x K matrix: 

{Zij i > j, 
i = j, 
-z^ i < j. 

The right hand side of expression fj4*3|) can be neatly written as a matrix integral: 

/ i \ N 
tt-^ J {K \{dZM)e-^ ZZl Pf( % C zX , (44) 
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where Q^ K > = {Z E Q KxK | Z T = —Z T } is the space of skew-symmetric complex 
matrices, X(Z, Z^) is the Lebesgue measure on Q^ K ' as described above. Note that 
the dimension of the integral in the right hand side of (1441) is iV-independent. The 
size of the original matrix only enters the integral as the power of the Pfaffian in 
the integrand. This allows one to calculate the large iV-limit of (jSJ) using the 
Laplace method. To facilitate the application of asymptotic methods, we re-scale 
the integration variables using (Z, Z^) — > yN(Z, Z^), which gives 

K(K-l) N K Wif K(K-l) 

= 7i s — 2~ — N—N — § — J N (45) 




N 

\{dZ,dtf)e-^ TrZZ 'Pf\ " V^ X 

The integrand in is now of the form exp(iV\F/v(Z)), where Fjy is a slow function 
of N. 

Applying the Laplace method. We will show that the integral localizes 
onto the subset 

CM = {Z E Q {K) | ZZ^ = I}. (46) 
To do this it is convenient to split J into two parts, 



Jn,o — / Jn,i 

JQWns JqW\s 

where 

S = {Z | fi k (Z) E [1/2, 2] for k = 1, . . . , K} 

and (fik{Z) : 1 < k < K} are the singular values of Z. Note that S is compact 
and contains only non-singular matrices. We first bound the integral Jjv,i, aiming 
to show that it is of smaller order that Jn,o- For c > write (Afc(c) : 1 < k < 2K) 
for the singular values of the matrix 

Z cX 
-cX 

One may bound the difference of two sets of singular values (/ik{A)) and (fik{B)) 
via \fji k (A) -Hk(B)\ < \\A-B\\. Thus 

lAfcfc) - A* (0)1 < \\cX\\ = ex* 
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where I. Then 

Pf 



Z cX 
-cX Zt 



det 



2K 



< 



Z cX 
-cX Zt 

1/2 

fc=i 

2if 

JJ (A fc (0) + cx*) 
k=i 

K 

Y\_(fj,k(Z) + cx*). 



1/2 



1/2 



Using this in Jjv,i we find 



<7zV,l < 



\(dZ,dtf)e— 



k=l 



K 



n (/i fc (z) + (2/iv) 1 / 2 ^) 



fc=i 



jv 



JQ( K )\S 



where 



H N (z) 



1 1 

2 ~ N 

1 -v2 



In (z + (2/N) 1 / 2 x,). 



Note that H N (z) — > #(2) = |-2 2 — ln(z) and that -ff(-z) has a minimal value of 
i?(l) = |. On \ S there must exist at least one singular value \x lying outside 
[|,2], and for this value > | + 25 for an easily calculated S > 0. Thus for 

large TV, when Z G \ 5, we have 



K 



K 



k=l 



and the value of J/v,i is bounded by C(K)e NK / 2 e NS . This is is exponentially 
smaller than that of Jno, which we will see is, to leading exponential order, 
0{e~ NK l 2 ). 

Next, we will calculate the asymptotic expansion of Jn,o for large N. The Nth 
power of the Pfaffian in the integrand of Jjv,o can be simplified using the Taylor 
expansion for the Pfaffian: 



Pf(A) ~ ' - l ' L ' 



2VN 



(47) 



+ ^ {TrBA' l TrBA- 1 - 2TrBA~ 1 BA~ 1 ) + O 
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While we cannot pinpoint the exact reference for the original derivation of the 
above expansion, it can be easily derived using the Berezin integral representation 
of the Pfaffian. It is interesting to note that unlike the analogous determinant 
expansion formula, the series (147p contains finitely many terms. Using (HTl) . and 
the fact that the terms with Tr{BA~ 1 ) are zero in our case, we can re-write the 
iVth power of the Pfaffian in the integrand of J^,o as follows: 

2 V \ / \ N 



I Z / — JC \ / i \ 

Pf N y_rj x V ^ f J = det(ZZt)* (l + ±Tr(ZlXZX)+0(N-*)) 

= det(ZZ^e Tr ( zixzx ) (1 + OiN- 1 )) 

This allows us to express J^.o i n a form well suited for the application of Laplace 
formula: 

J N . = [ A(^,^t)e-^( T ^ t - lndct ( zzt ))e T K^ t ™)(l + 0(iV- 1 )). (48) 

The fact that S does not contain degenerate matrices, and the compactness of S 
allows one to pass the correction term 0(N~ l ) through the integral. In the limit 
N — > oo, the main contribution to ( 148]) comes from the neighborhood of the points 
of global minimum of the function 

K 

F{Z) = TrZZ ] - In det (ZZ ] ) = ^ (l4(Z) - 2 hx(p k {Z))) . 

k=l 

The global minimum value of F is K and it is attained on the set C^ K ' in (f46|) . 
namely the skew-symmetric unitary K x K matrices, which is a smooth sub- 
manifold of Q {K \ We will show that is a non-degenerate critical set, which 
means that the Hessian of F has the maximal possible rank at every point of C^ K \ 
Therefore we can use Laplace theorem [7] to calculate the asymptotic expansion 
of Jn,o- let (w, y) be local co-ordinates on Q( K ) such that the sub-manifold C^ K ' 
is locally determined by the set of equations y = 0; then 

f A(dZ, dZ^eS (TrZZt-lndct(ZZl)) e Tr( Z iXZX) 

JQ( K )ns 

/ , \ dim(Q( K ))-dim(c( K >) r 



e -NF\ 



Here fi(dw) is the measure on C^ K ' generated by the embedding C^ K ' C Q^ K ' and 
integration over transverse co-ordinates y. Explicitly, 

_ . *m(OW) 

»m= i pKy ~ uj n d Wk , (so) 

'det ffess(F) |g° (tw) *=i 
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where p(w,y) is the density of Lebesgue measure A(dZ,dZ^) with respect to local 
co-ordinates (w,y), the Hessian is defined as the matrix of second derivatives with 
respect to transverse co-ordinates y. In writing ()49l) we used the fact that the 
critical manifold lies a positive distance away from the boundary of fl S. 
Noting that F takes the value K on C^ K > and that 

dim(Q w ) - dim(C w ) = K{K - 1) - (dim(U(K)) - dim{Sp{K))) 

= K(K -1) - K 2 + ^K(K +1) 

= \k{k-i), 



we reach 



J N0 = e- E f(2rcN)-^ / p{dw)e T < z ' xzx ) (1 + O(N^)). (51) 

Jew 



Collecting together (J42]), (g5]) and ()5l]) we find 

K r 

x 2 ,..., x K ) = c 2 (N, K)A(x) TT e~ x * \ 

where 



p N (xi, x 2 ,..., x K ) = c 2 (N, K) A(x) I I e""* / //(e^e 1 ^™) (1 + (1)) 



n/ N-k\ K(K-l) (N-K)K 
(\S N ^\7T-—) 

fc=l 



2 



(N - K) — 2 — (AT - K) — a — e 2 — (2tt(A^ - AT)) ~ 

and C(K) denotes a constant only depending on K. It is lengthy but straightfor- 
ward to check that c 2 {N,K) — > c%(K) > as A" — > 00 and hence that limiting 
modified density p(xi, x 2 , ■ ■ ■ , xk) = liniAr_>oo p~n(xi, x 2 , . . . , xk) exists and is given 
by 



K 2 f 

p(x u x 2l ...,x K ) = c 3 (A:)A(x) TT e -x * / 



M^)e Tr ( ZtXZX ). (52) 

The explicit value of c^{K) will be determined later by using properties of the 
densities p. In the next subsection we will find a parameterisation of the integral 
in the right hand side of (|52|) . which will allow us to calculate it very efficiently 
using the standard tools of random matrix theory. 



Recasting as an integral over the unitary group. An important property of 
the function F is its invariance with respect to the following action of the unitary 
group U{K) on Q^: 

U{K) x U(K) 

(U,Z) i-5> UZU T eU(K). (53) 
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Namely, for any A e U(K), we have F(AZA T ) = F(Z). The decomposition 
theorem for skew symmetric unitary matrices [H] states that 

Z = UJU T } (54) 

where U is a unitary matrix, J is the canonical symplectic matrix. Notice that 
( 1541) does not determine the unitary matrix U uniquely: indeed Z—tZifU US, 
where S is a unitary matrix satisfying SJS T = J. The set of such matrices is a 
subgroup of U(K) called the symplectic group Sp(K), that is 

Sp(K) = {S e U{K) | SJS T = J}. 

It can be checked that the critical manifold C^ K ' can be identified with the factor 
space of U(K) with respect to the action of Sp(K) on U (K) via right multiplica- 
tion, that is 

C (K) ^ u( K )/Sp(K) (55) 

The [/(fO-action ([53]) on QW preserves the critical manifold and induces the 
[/(iC)-action on C. Using the parameterisation (154]) of C^ K ' this induced action 
can be written explicitly: 

U(K) x C {K) -> C {K \ 

(A,[U\) ^ [AU], (56) 

where [U] is an equivalence class of U G U(K) with respect to right multiplications 
by elements of Sp(K) C U(K). In the vicinity of a critical point Z c e C^', 

F(Z c + ^) = / f+ irr(^t +W+ ... (57) 

We notice that the quadratic form describing the second order term in the above 
Taylor expansion of F is U(K) -invariant and has the maximal possible rank equal 
to \K{K — 1) = dim(Q^- ) ) — dim(C( x )). We rewrite the integral from ( 152]) . using 
the mapping (155]) . as 

/ ^(dw)e Tr ( ztxzx ) = [ ft(dU)e Tr ( z ^ xz ^ x )) (58) 

Jem Ju(K)/Sp{K) 

where Z(U) is given by ( 1541) and £i(dU) is the pull back of the measure // on 
the critical manifold. We can work out an explicit expression for /i in local co- 
ordinates on using the general formula (150]) . We will not do this; instead we 
will characterise p, up to a multiplicative constant by establishing its symmetry 
with respect to the ?7(iV)-action on C^ K \ Recall that the measure /j, is determined 
by the Lebesgue measure on Q( K > and the determinant of the quadratic form in 
the right hand side of (l57j) . It is easy to check 
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(i.) The Lebesgue measure A and the quadratic form Tr(5ZZ^ + ZSZ^) 2 on 
are invariant with respect to the £7(-fT)-action (|53|) . 

(ii.) The critical manifold C^ K > is invariant with respect to the U "(i^-action. 

(iii.) The restriction of the quadratic form Tr(5ZZ^ + ZSZ^) 2 on to 
has maximal rank. 



A calculation employing elementary tools of differential geometry [5] shows that 
the above three observations imply the invariance of the measure \x with respect 
to the induced action of U(K) on the critical manifold U(K) / Sp(K) defined by 
( 1B"6l) . Therefore fi is a Haar measure on the symmetric space U(K)/ Sp(K), which 
is unique up to normalization. It is generally easier to work with integrals over the 
whole unitary group rather than a factor space. As we have established already, 
the measure ft is invariant with respect to the action of U(K) on C^ K '. Note that 
the function Tr (Z(UYXZ(U)X^ which determines the integrand of fl58|) is also 
[/(-fT)-invariant. Therefore, by Weyl's theorem [10] Chapter X, 

[ fi(dU)e T < z ^ xz ^ x ) = [ f, H (dU)e Tr ( z ^ xz ^ x ) (59) 

JU(K)/Sp(K) Ju(K) 

where /i# is an appropriately normalized Haar measure on the unitary group. 
We will determine the normalization factor later using the properties of spin-spin 
correlation functions. Substituting ( 1541) into ( 1591) we find that 



[ fx H (dU)e T < z ^ xz ^ x ) = [ MdU)e- Tr ( JHJHT ) (60) 

JU(K) Ju(K) 

where if is a Hermitian matrix with eigenvalues Xi,X2, ■ ■ ■ ,Xk given by H = 
UXU^ . Tracing back we find that the large N behaviour of 

K 

E Y[ det (Mi - x m ) 



m=l 



- the expected value of the product of K characteristic polynomials in the real 
Ginibre ensemble - turns out to be determined by the integration of the symplectic- 
invariant Gaussian weight exp(—TrJHJH T ) with respect to unitary degrees of 
freedom. See [18] for a discussion of the origin of the connection between real 
Ginibre and symplectic ensembles. 

We may rewrite the integral in (|60|) as 

MdU)e T < HHR ) (61) 



U(K) 

where H = UXU^ and 

H R = JH T J T 



27 



is a 'symplectic' involution on the space of complex K x K matrices, see [15] for 
details. Following Mehta, we will call matrix M self-dual if M = M R and anti- 
self-dual if M = —M R . It is easy to check that any even-dimensional matrix can 
be uniquely represented as a sum of a self-dual and anti-self-dual matrices. Let 
ASD {K) be the linear space of all anti-self dual K x K matrices. In order to 
perform the integration over the unitary group in the right hand side of (l6Tj) . we 
use use the following transformation found in |17] . |16j : 

eTrHH n = f A ( dA ) e TrAi+2V2TrHA^ (g2) 

Jasd( k ) 

where A(gL4) is a Lebesgue measure on Zk is a normalization constant. 

The absolute convergence of the above integral can be checked using the decom- 
position theorem for anti-self-dual matrices, see [H]: for any A e ASD^ there 
exists V G U(K) and a diagonal matrix G with entries ±#i, ±# 2 , ■ ■ ■ i 9k , where 
9i > 0, 2 > 0, . . . e K/2 > 0, so that 

a = ivevK 

Substituting fl62l) into (!6T|) and using the invariance property of the Haar measure 
Hh we get: 

f m {dU)e Tr ( HHR )=C K e^*l f ^ H {dU) [ v(d6)e- Tr ( @2 ) e t2V ~ 2TrUXUt& 

JU{K) JU(K) JR^ 2 

where v(d9) is the measure on the eigenvalues of unitary matrices induced by the 
mareinalization of the Lebesgue measure on over the unitary degrees of 

freedom, explicitly [15J 

Kjl 

v{d9) = a (±o u ±e 2 , . . . ± 9k ) Y[ e k de k . 

k=l 

We now allow the constants Ck, depending only on K, to change value from line 
to line. We can integrate over the unitary group U(K) using Harish-Chandra- 
Itzykson-Zuber formula [9], [13]. The result is 

I JTT\ Tr(HH R ) 

U{K) 

r K/2 

= C K e^ x l Aix)- 1 / n 6ke ~ el d9k det e i2Vlxje " -l<j,l<K 

The remaining integration over the singular values 61,62, ■ ■ ■ ,&k/2 is carried out 
using the de Bruijn formula [1], yielding 



/ VH(dU)e Tr ( HHR )=C K e^A(xrPf 

'U(K) 



yXi XjJC 



l<i,j<K 
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Combined with ( 152|) this yields 



p{xi,x 2 , . . .,x K ) = C k Pf 



yXi XjJC 



l<ij<K 



(63) 



The correlation functions of spins can be computed by integrating p N with respect 
to space variables: 



P~N(yi,y2, 



Vk) 



(64) 



This leads to the spin-spin correlation 

A' 



El[s kk (M l ) = C k Pf 



k=l 



r 2z dz 



(65) 



l<i,j<K 



The constants C k can be found inductively in k by allowing x 2 fc I %2k-i, and noting 
that p(x 1 ,Xi) = 1. This yields Ck = (8/ir) K /' i . Expression ( 165|) coincides with 
the single time spin-spin correlation for a system of one dimensional annihilating 
Brownian motions under the maximal entrance law, see [22] for details. As shown 
in [22], substituting fl65|) into formula 02]) to obtain we get the first statement of 
Corollary 9 of 0. 
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